REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply 
with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS.  _ ■ _ 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

1 1  -06-2003  Final  Report  27  July  2002  -  27-Jan-03 


4.  TITLE  AND  SUBTITLE 


Accurate  Stress  Intensity  Factor  Solutions  for  Unsymmetric  Corner  Cracks  at  a 
Hole  Subject  to  General  Loading 


5a.  CONTRACT  NUMBER 

F61775-02-WE023 


5b.  GRANT  NUMBER 


5C.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


5d.  PROJECT  NUMBER 


Dr.  Borje  Andersson 


5d.  TASK  NUMBER 


5e.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

The  Swedish  Defense  Research  Agency 

Ranhammarsvagen 

Stockholm  SE-172  90 

Sweden 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


EOARD 

PSC  802  BOX  14 
FPO  09499-0014 


11.  SPONSOR/MONITOR’S  REPORT  NUMBER(S) 

SPC  02-4023 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


20040715  184 


14.  ABSTRACT 


This  report  results  from  a  contract  tasking  The  Swedish  Defense  Research  Agency  as  follows:  A  basic  research  program  will  be  conducted  to 
calculate  stress  intensity  factor  (K)  solutions  not  only  for  symmetric  but  also  unsymmetric  corner  cracks  at  a  hole  subject  to  general  loading. 
The  p-version  finite  element  method  (FEM)  with  a  mathematical  splitting  scheme  will  be  used  to  enable  efficient  and  accurate  calculations.  All 
structurally  significant  crack  shapes  are  to  be  considered.  In  addition,  combinations  of  crack  depth  to  crack  length  (a/c),  crack  depth  to  sheet 
thickness  (a/t),  and  hole  radius  to  sheet  thickness  (r/t)  ratios  are  to  be  analyzed  at  each  side  of  the  hole;  thus  more  than  106  solutions  are  to 
be  developed  with  control  of  the  error  in  the  computed  K  solutions.  The  loading  conditions  will  be  remote  tension,  remote  bending,  and  pin 
loading  (bearing).  The  new  K  solutions  will  be  implemented  in  AFGROW. 


15.  SUBJECT  TERMS 

EOARD,  Materials,  Fracture  mechanics,  Failure  Mechanisms,  Metallic  Materials 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

UNCLAS 

b.  ABSTRACT 

UNCLAS 

c.  THIS  PAGE 

UNCLAS 

17.  LIMITATION  OF 
ABSTRACT 

UL 


18,  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 
OF  PAGES  CHARLES  H.  WARD,  Lt  Col,  USAF 
4 

19b.  TELEPHONE  NUMBER  (Include  area  code) 
+44  (0)20  7514  3154 


Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39-18 


17. FEB. 2004  17=01 


FOI  BROMMfl  46  8  2534S1 


NK.  /7k3 


'a.d./Z 


(1)  In  accordance  with  Defence  Federal  A<  :quisition  Regulation  252.227*7036,  j 

Declaration  of  Technical  Data  Conformity  (Jan  1997),  ! 

“The  Contractor,  Swedish  Defense  Research  Agency,  hereby  deleaves  that,  to  the  best  of 
its  knowledge  and  belief,  the  technical  data  delivered  herewith  under  Contract  No. 

F6 1775-02- WE023  is  complete,  accurate,  anji  complies  with  all  requirements  of  the 
contract. 

DATE:  1  f-WU , 

Name  and  Title  of  Authorized  Official: 

1 

Ingvar  Roos 

(2)  In  accordance  with  the  requirements  i  i  Federal  Aquisition  Regulation 
52.227-13,  Patent  Rights — Acquisition  by  i  he  V.S.  Government  (Jun  1989), 

(B)  “I  certify  that  there  were  no  subject  inventions  to  declare  as  defined  in  FAR 
52.227-13,  during  the  performance  of  tlqs  contract,” 

DATE:  j  ,  2a 

Name  and 


1  i-weu  ,  _ 

Title  of  Authorized  Official-^: 


Ingvar  Roos 


DTIC  Copy 


/}Q  Fctt-cft-  /oo4 


SUMMARY  REPORT 


Accurate  Stress  Intensity  Factor  Solutions  for  Unsymmetric  Corner 
Cracks  at  a  Hole  Subject  to  General  Loading 

EOARD  Supported  Research 

Reference:  SPC  024023 


Date:  Feb  26, 2004 


Background 

In  fatigue  design  of  aircraft  structures  so  called  stress  intensity  factor  data  (K)  are  needed 
for  cracks  of  various  sizes.  One  of  the  most  important  geometries  is  small-cracks  located 
at  rivet  holes  in  the  aircraft  skin. 

It  was  recently  shown  that  the  few  (i.e.  of  order  hundred)  existing  numerical  solutions 
have  errors  ranging  from  2%  up  to  40%,  which  is  unacceptably  large  [1]. 

The  need  for  accurate  K- data,  and,  for  a  larger  parameter  range  motivated  the  present 
research. 

In  the  present  work  (Ref:  SPC  024023)  about  5  million  new  solutions  have  been  derived 
with  a  guaranteed  relative  error  below  1%.  Novel  mathematical/numerical  methods  were 
used  in  the  project. 

Crack  geometries  considered 

The  figure  below  shows  the  geometries  studied  (Note:  wedge  angle  used  is  0  degrees 
and  b=t).  There  are  two  quarter  elliptical  cracks  at  the  hole  boundary.  More  detailed 
technical  results  from  the  present  project  are  reported  in  [2].  For  a  method  description, 
see  [3]. 
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Figure  1  Two  Unsymmetric  Corner  Cracks  at  a  Cylindrical  Hole 


Hole  Geometry 

Cylindrical  hole  with  radius  R=D/2  (countersunk  angle  =  0°) 
b/t  =1.0 


Part-Through  Elliptical  Crack  Sizes 

The  following  parameter  sets  were  analyzed. 

a/cj  =  0.1, 0.111, 0.125, 0.1428,  0.1667,  0.2,  0.25,  0.333,  0.5,  0.667, 0.75,  0.80, 1.0, 1.25, 
1.333,  1.5,  2.0,  3.0, 4.0,  5.0,  6.0,  7.0,  8.0,  9.0,  10.0 

a/t  =  0.1, 0.2, 0.3, 0.4, 0.5,  0.6,  0.7,  0.8,  0.9,  0.95,  0.99 

r/t  =  0.1, 0.111,  0.125, 0.1428,  0.1667, 0.2, 0.25, 0.333,  0.5,  0.667,  0.75, 1.0, 1.25,1.333, 
1.5, 1.666, 2.0,  3.0, 4.0,  5.0,  6.0, 7.0,  8.0,  9.0, 10.0 

In  the  present  work,  R=2  units  and  W=200  units  were  used  in  the  numerical  calculations. 
For  other  dimensions,  all  AT-data  scales  as  the  square  root  of  (r/2),  where  r  is  the  actual 
radius. 

Material  data  in  the  plate  is  isotropic  with  Poisson’s  ratio  0.3. 

For  25  r/t  ratio’s,  1 1  a/t  ratio’s,  25  different  a/c  ratio’s  and  three  types  of  loading  (tensile, 
bending  and  wedge  loading),  the  total  number  of  K  -solutions,  in  case  of  two  cracks,  is 
25*11*25*25*3=  5156250.  The  number  of  single  crack  solutions  is  20625.  For  each  of 
the  over  5  million  solutions,  the  mode  I,  mode  II  and  the  mode  III  stress  intensity  factors 
have  been  computed  at  77  points  at  each  crack  front  (hence,  over  2  billion  AT-numbers  are 
stored  in  the  stress  intensity  factor  library). 


Stress  Intensity  Factor  Library 


All  computed  X-data  are  today  stored  on  the  computer  strato2  at  the  USAF  Academy  in 
Colorado  Springs,  CO. 

Summary 

The  work  stipulated  in  the  contract,  Order  Number  F61 775-02 -WE023,  has  therefore 
been  completed. 


Acknowledgements 

A  part  of  all  computations  (about  40%,  corresponding  to  about  200,000  CPu-hours)  was 
supported  by  a  grant  of  computer  time  from  the  DOD  High  Performance  Computing 
Modernization  Program  at  ERDC. 


Sincerely, 

Borje  Andersson 

Research  Director,  Structures&Materials 
The  Swedish  Defence  Research  Agency 
SE  172  90  Stockholm 
Sweden 

References: 

[1]  s  Fawaz,  B  Andersson,  ’’Accurate  Stress  Intensity  Factor  Solutions  for  Unsymmetric 
Comer  Cracks  at  a  Hole”.  Presented  at  Fourth  Joint  DoD/FAA/NASA  Conference  on 
Aging  Aircraft,  May  15-18, 2000,  MO  US 

[2]  S  Fawaz,  B  Andersson,  and  J  C  Newman  Jr,  ’’Experimental  Verification  of  Stress 
Intensity  Factor  Solutions  for  Corner  Cracks  at  a  Hole  Subject  General  Loading”, 
International  Committee  for  Aeronautical  Fatigue,  ICAF  2003,  Lucem,  Switzerland,  May 
5-7, 2003 

[3]  S  Fawaz,  B  Andersson,  ’’Accurate  Stress  Intensity  Factor  Solutions  for  Comer  Cracks 
at  a  Hole”,  Engineering  Fracture  Mechanics,  Vol.  71/9-10,  pp  1235-1254, 2004 


